High accuracy measure of atomic polarizability in an optical lattice clock 
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Despite being a canonical example of quantum mechanical perturbation theory, as well as one of 
the earliest observed spectroscopic shifts, the Stark effect contributes the largest source of uncer- 
tainty in a modern optical atomic clock through blackbody radiation. By employing an ultracold, 
trapped atomic ensemble and high stability optical clock, we characterize the quadratic Stark ef- 
fect with unprecedented precision. We report the ytterbium optical clock's sensitivity to electric 
fields (such as blackbody radiation) as the differential static polarizability of the ground and excited 
clock levels adock = 36.2612(7) kHz (kV/cm)"'^. The clock's fractional uncertainty due to room 
temperature blackbody radiation is reduced an order of magnitude to 3 x 10^^^. 

PACS numbers: 32.10.Dk, 32.60.+i, 06.20.fb,44.40.-|-a 



An atom immersed in an electric field i?a becomes 
polarized — the electronic wave-function is stretched into 
alignment with the field. Generally, energies of the 
lowest- lying electronic quantum states \i) are reduced 

by —^^ctf^ El (see Fig. [ij) where a-°^ is termed a state's 
static polarizability [T]. The scaling of the Stark effect 
is quadratic because i?a is responsible for inducing, and 
also interacting with, an atomic dipole moment. 

Neutral atom lattice clocks [5] employ lO'^-lO^ ultra- 
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FIG. 1. a) A scale drawing of transparent conductive elec- 
trodes surrounding a one-dimensional optical lattice of ytter- 
bium atoms. Four pairs of metallic pads (one of which is 
highlighted) allow m situ interferometric measurement of the 
electrode spacing d. b) Deviations cr_B of the applied electric 
field Ea, from the ideal field created by infinite-planar elec- 
trodes become negligible in the central 1 cm region (dashed 
white circle), c) Relevant energy levels, laser transition wave- 
lengths, and linewidths (r/27r) of ^^^Yb. Clock states ^5*0 and 
^Po are shown Stark-shifted by Ea- 



cold alkaline-earth atoms tightly confined in an opti- 
cal standing wave potential so their intrinsically nar- 
row, largely imperturbable, ^Sq ^ ^Pq transitions [3] 
may establish stable and accurate frequency and time 
references [H [5] . In analogy to a pendulum's oscillation 
slowing due to thermal expansion, a ytterbium lattice 
clock slows when electrically stretched, or polarized, by 
thermal blackbody radiation (BBR) fields [7]. This 
phenomenon has been measured in the cesium foun- 
tain primary standard [8] ^ and other optical transi- 
tions [3 [H EI] . 

No shield at finite temperature protects a clock atom 
from the time varying electric field of BBR, the electro- 
magnetic energy absorbed and re-emitted by all matter in 
thermal equilibrium according to the Stefan-Boltzmann 
law [12] . Inside a hollow shell of opaque matter (a black- 
body), the time-averaged electric field intensity depends 
strongly on temperature: 
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where a « 1/137 is the fine-structure constant, fee is 
Boltzmann's constant, and T is the blackbody's absolute 
temperature [13] . Near room temperature, the spectrum 
of radiation is peaked strongly near 9.6 /im — invisible 
to the eye and also far detuned from strong electronic 
transitions in ytterbium. 

The ytterbium clock frequency {v 518 THz) is 
shifted by the net BBR Stark effect of the two clock 
states, which can be expressed as 



AfBBR 



where a'g}e are the static polarizabilities of ground and ex- 
cited states (-^^o and '^Pq, respectively), and ?/ciock(P) C] 
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is a small computed parameter accounting for the dy- 
namic aspect of the BBR field [Mj. At room tempera- 
ture, r?ciock(300 K) = 0.0145(15) is a small contribution 
to Ai^BBR ~ —1-3 Hz. More significantly, knowledge of 
aciock = ctf^ — c^g°^ is theoretical, and limited to 10% 
accuracy due to the complexity of this many-electron 
atom [Hlllllle]. 

To measure aciock: and reduce the clock's uncertainty 
due to BBR, we fitted electrodes [17] to an existing 
ytterbium clock apparatus [B]. A voltage V on ideal 
electrodes spaced by d in vacuum creates an electric 
field Eg, = V/d, shifting the clock transition by Ai/ = 
— ^ciciockiy / d)^ ■ Deviations from this infinite-parallel 
plane capacitor model are bounded at the 10~^ (1 ppm) 
level by designing a large electrode diameter-to-spacing 
ratio, ensuring a high degree of parallelism, and centering 
the atoms radially within the electrodes. Perturbations 
due to dielectric and conducting mounting structure con- 
tribute similar amounts of field uncertainty. 

The electrodes, shown in Fig. [l] are comprised of 
rigidly spaced parallel fused silica plates (101.6(1) mm 
in diameter) featuring a transparent conductive coat- 
ing on the inner surfaces. The outer surfaces are anti- 
reflection coated for all relevant laser wavelengths. The 
electrode separation d = 15.03686(8) mm is maintained 
by three fused silica rods bonded 45(1) mm from the 
center axis with hydroxide catalysis [HI [19]. d is de- 
termined interferometrically by measuring {in situ) the 
free-spectral-range t'fsr = c/2d of planar etalons formed 
by 90%-reflective metallic pads (33 nm gold on 2 nm 
chromium) deposited over the 0.3 nm indium-tin-oxide 
inner electrode faces [17]. Each 6 mm square pad is off- 
set 28 mm from center. An external cavity diode laser 
is tuned by 17 THz around Ap — 766 nm to observe 
~ 10 of N{ ^ 1700 etalon transmission peaks. Each 
transmission feature, located with a wavelength meter 
to ±50 MHz, has a linewidth of 500 MHz, consistent 
with a finesse « 20. Spacing the observed peaks loga- 
rithmically allows an efficient least-squares determination 
of Ufsr- Systematic wavelength and fringe-center inaccu- 
racy is largely divided down by iVf. Uncertainty in the 
metal pad thickness contributes 3x 10^^ uncertainty to d. 
Gold's index of refraction (n^ « 0.16) [20] varies gently 
about Ap; variations in visr with fringe index due to mir- 
ror phase shifts contribute an error of 1.2 x 10^^. Stray 
etalons add a similarly sized line-pulling error. Electrode 
parallelism is constrained by the measured finesse as well 
as spatially independent pad-pair measurements of d. 

Thin strips of silver-loaded epoxy join insulated wires 
to electrode perimeters. Two wires are redundantly 
bonded to each electrode to establish the sheet resistance 
of the conductive layer (J?ito = 3 kil) . The parasitic resis- 
tance between the electrodes and the grounded vacuum 
structure is i?icak = 316(9) (i?icak = 14.3(5) Gn) 
below (above) an observed field-emission threshold oc- 
curring near 800 V. Leads to the electrodes have current- 
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FIG. 2. Integrated clock laser error signals are shown for 
iSa ~ 700.6 V/cm under three interleaved conditions: both 
electrodes grounded ('A', grey dots), and each electrode at 
high voltage with the other grounded ('B' and 'C, blue and 
red dots). Laser drift, common to these signals, is removed 
from the quadratic Stark shift signal (black open circles, see 
Eq.|3| and the stray field signal Aub^Auc (blue open circles). 
Dashed lines show standard deviations for these signals. A 
passive feed-forward linear drift canceler reduced clock laser 
drift to 4 Hz over this 1000 s data run. 



limiting resistances of 10 kfl each; the worst-case i?a er- 
ror from voltage division is 1.4 x 10^^. We constructed 
a regulated voltage source producing 100 V-1050 V with 
1.0 X 10"^ instability over 1 s-1000 s EH- 

A clock interrogation cycle (360 ms) begins with slow- 
ing, cooling, and trapping ^''^Yb from an atomic beam 
with a magneto-optical trap operating first on the ^Sq -(->■ 
^Pi transition (399 nm), then on the narrower ^5*0 O ^Pi 
transition (556 nm). Ultracold atoms (10 /iK) are con- 
fined by an optical lattice (ID) at the so-called magic 
wavelength near 759 nm resulting in no net ac-Stark shift 
between ^5*0 and ^Pq [HI HI] . Atoms are optically pumped 
into one nuclear-spin state (tof = ±1/2) of ^Sq. A tt- 
polarized pulse (100 ms) of resonant 578 nm laser radia- 
tion coherently excites one of two Amp = transitions 
(split by 500 Hz by an applied magnetic field B) to the 
long-lived ^Pq state (see Fig. [l|. A series of laser pulses 
converts the resulting clock state populations into fluo- 
rescence signals which are then normalized against atom 
number fluctuations [6]. 

The stabilized 578 nm laser is independently locked 
to the atomic transition under three interleaved con- 
ditions: both electrodes grounded (condition 'A'), and 
each electrode at high voltage with the other grounded 
('B' and 'C'). Each condition, mediated by opto-coupled 
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reed-relays, lasts Tv = 2.9 s, during which the clock laser 
maintains a fractional frequency stability approaching 
3 X 10~^^ [23]. In each period Tv, two interrogations are 
performed on each side of both nuclear-spin spectroscopic 
features. Slow laser drifts are common to all three inte- 
grated error signals, Ai/^, Avb, and Ai^c (see Fig. [2]). 
The quadratic Stark shift is 
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Aiy = i(Ai/B -I- Ave) - I^VA- 



(3) 



The total Allan deviation [24J is used to determine the 
statistical uncertainty of Av. For presented data, || 
-^lattice (both were perpendicular to i?), though other 
configurations were examined. 

Reversing |17j yields information about stray elec- 
tric fields E's parallel to or differential contact po- 
tentials ^25] (e.g. one electrode may develop a thin layer 
of ytterbium deposition). The difference Avb — Ave — 
2aciockEa- Es reveals that \Es-z\ « O.I V/cm (see Fig.|3|; 
temporal drift and weak correlation with _Ea are observed. 

Shifts due to a truly static field Eg subtract completely 
in Eq. |3] However, time dependent changes, notably 
those correlated with the polarity of E^, do not. Be- 
cause of careful experimental design, we are not aware 
of any appreciable stray field source with such a corre- 
lation. Nevertheless, because an increase in dwell time 
Tv potentially allows an accumulation of unknown stray 
charge (and thus a correlated stray field), we varied Tv 
over 0.8 s-I2 s and resolved a small but negligible cor- 
relation in Av at applied fields twice the maximum used 
for reported data. We observed no systematic variation 
in the measured polarizability adock with E^^. We note 
that for data presented here, i?a remained three orders 
of magnitude below the dielectric strength of fused sil- 
ica, and five orders of magnitude below the characteris- 
tic level [26] for ITO electron emission. The time con- 
stant for electrode charging is 8 ^s. Typically, 100 ms 
is allowed for settling. Connecting high voltage to both 
electrodes creates an electric-gradient field; by observ- 
ing Av < 30 mHz with 2 kV applied, we constrain the 
atoms' radial position to ±10 mm, consistent with visual 
observations. 

Fig. [3^ shows the observed clock frequency shift 
quadratically as a function of E^,. When fit to a polyno- 
mial, the data are consistent with no quartic, cubic, lin- 
ear, or offset terms — an ideal demonstration of the Stark 
effect as non-degenerate perturbation theory. No inho- 
mogeneous line broadening is observed with increased 
shift, so the fractional statistical uncertainty in Av re- 
duces as E^^. In contrast, the uncertainty of the applied 
voltage (the dominant systematic uncertainty) rises as 
E^ according the specifications of our commercial volt- 
meters. Fig. [3]; plots the polarizability inferred at each 
i?a- Table |T] lists the sources of measurement uncertainty 
at a particular applied field. Taking the mean of all mea- 
surements, weighted by the total standard errors, we de- 
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FIG. 3. The clock slows when stretched, a) Residuals of a 
quadratic fit Av = — |Qciocki5f display measurement uncer- 
tainties in Av (below) and Ea, (right), b) An inset depicts 
lattice trapped atoms and the relative orientations of -Eiattice, 
Ea, and hypothetical static stray field Eg. c) At each E^, 
we display Qdock = —2Av/E^ with combined statistic and 
systematic uncertainties (see Table [l|. An inset shows data 
at higher resolution. Solid and dashed lines show the final 
result and standard error, respectively, d) The component of 
a stray field Eg \\ E^ is precisely detected upon E^ reversal. 
We observe more temporal variation in Eg than correlation 
with Eb,. Data with different marker styles were acquired on 
separate days. 



termine adock = 36.2612(7) kHz(kV/cm)-2. A least- 
squares functional fit (Fig. |3^) yields a consistent value 
for Qfciock- Table [H] demonstrates the agreement between 
this measurement and theoretical predictions. 

Neither static nor BBR fields cause spin-magnitude- 
dependent (cx |mi?p) tensor Stark shifts because both 
clock states have insufficient total angular momentum 
(F = F' = 1/2) [27]. Spin-sign-dependent (cx mp) vec- 
tor Stark shifts are absent as well: BBR has no net po- 
larization and static fields lack the time dependence to 
be circularly polarized [2S]. No opposite-parity states 
lie close to either clock state so no linear dependence 
of Av on i?a is expected or observed. A third-order ef- 
fect (29j mixing the polarization due to the optical lattice 
(-^'lattice ~ 10 kV/cm) and E^, is expected to cause a I0~^ 
fractional error at the highest i?a- A fourth-order term 
Av cx E^ (the hyper-polarizability) is responsible for a 
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TABLE I. Uncertainty budget for a representative datum. 
Errors in Ea, contribute twice the uncertainty as those in Av 
due to the dependence adock = —^Aiz/E^. This factor of 
two is included in the tabulated quantities below. The total 
uncertainty in ciciock, for the particular measurement shown, 
is found by summing all contributions in quadrature. 



Uncertainty source 


xlO~^ 


Notes /conditions 


Shift statistical error 


8.3 


Ai^ ^ -3603.77(3) Hz (1800 s av- 
eraging) 


Higher-order Stark shifts 


0.01 




Electric field (E^) errors: 




-C/a — 44o.ood(_4_) V/cm 


Voltmeter systematic 


16.4 


Regulated 670.3966(55) V 


^leak voltage division 


0.1 


-iicak — ^-J- nA; zU KSi Icaas 


Finite electrode size 


1 


Atoms centered ±10 mm 


Electrode parallelism 


4 


^wcdge < 7 /irad 


Electrode deformation 


0.8 


Warping of fused silica by gravity 


Dielectric spacers 


2 


Perturbation of ideal field due to 
three fused silica posts 


Spacing d (statistical) 


1.6 


A'^f > 1700 fringes spanned 


Spacing d (systematic) 


9 


Fringe centering, wavemeter ac- 
curacy, stray etalons, stability 


Etalon probe tilt, </> 


0.3 


(1 — cos 0) error, (f) <C 0.5 mrad; 
retro-coupling single-mode fiber 


Yb thermal beam 


0.06 


Dielectric (e^ - 1) 8 x 10~^ 


Stray fields, static 


0.04 


Uncertainty in Eg^ reversal 


Stray fields, varying 


2 


correlation with 


Uncertainty in adock 


21 





ports with less than unity emissivity and opacity. An 
effective 1 K uncertainty in T leads to a fractional un- 
certainty of 3 X 10^^^ in the room temperature clock. A 
cryogenically shielded environment at T = 77(1) K with 
carefully controlled optical access can reduce the BBR 
shift uncertainty to the 1 x 10^^^ regime [3T]. A stray 
static field of 0.1 V/cm shifts the ytterbium clock tran- 
sition -0.18 mHz, a fractional change of 4 x 10^^^. A 
conductive enclosure at any temperature further ensures 
these, or smaller, stray static fields [52] . 

We note that our present uncertainty in adock is com- 
petitive with the best known atomic or molecular polar- 
izability, that of helium [33] . Beyond timekeeping, possi- 
ble metrological applications of the present work include 
high voltage measurement without the use of resistive di- 
viders, and more directly, an atomic electric field meter 
sensitive at moderate field strengths. 
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tributing to electrode fabrication, T. Fortier and S. Did- 
dams for femtosecond optical frequency comb measure- 
ments, and E.N. Fortson for a useful discussion on po- 
tential systematic effects. 



TABLE II. Comparison with theoretical predictions. Results 
are also presented in SI and frequently used atomic units 

(a-u.) m. 

— (0) (0) 



[kHz (kV/cm)"^] 


[ a.u. ] 


10-^3 [C mVV] 


Reference 


40.1(3.7) 


161(15) 


2.65(25) 




38.6(4.0) 


155(16) 


2.56(26) 


m 


33(13) 


134(51) 


2.21(84) 


pi 


36.2612(7) 


145.726(3) 


2.40269(5) 


this work 



similarly sized effect. We observed no systematic effect 
in Ai/ upon varying the lattice intensity or polarization. 
We ensured that the atomic density did not systemati- 
cally vary with application of E^] such a correlation could 
introduce contamination from the cold collision shift ^30] . 
We resolved no shift by systematically varying the elec- 
trodes between grounded and floating conflgurations. 

In an ideal blackbody environment at 300 K, the clock 
transition frequency is reduced from its value at T = K 
by Ai^BBR — —1.273(1) Hz, where the uncertainty is 
now dominated by that of jyciock (see Eq. [2|. In prac- 
tice, precise knowledge of adock and ?7ciock is not suffi- 
cient to determine the effect of BBR, since knowledge of 
the environment is limited. Non-uniformities arise due 
to temperature gradients, a hot effusive oven tip (850 K) 
and heated viewport (600 K), each about 30 cm away 
from the trapped atoms, and vacuum walls and view- 
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